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Abstract 

We further examine the capabihty of the 7 and 8 TeV LHC to explore the param- 
eter space of the p(henomenological)MSSM with neutrahno LSPs. Here we present 
an updated study employing all of the relevant ATLAS SUSY analyses, as well as 
all relevant LHC non-MET searches, whose data were publically available as of mid- 
September 2012. We find that roughly 1/3 of our pMSSM model points are excluded at 
present with an important role being played by both the heavy fiavor and multi-lepton 
searches, as well as those for heavy stable charged particles. Nonetheless, we find that 
light gluinos, lst/2nd generation squarks, and stop/sbottoms 400 — 700 GeV), as 
well as models with 1% fine-tuning or better, are still viable in the pMSSM. In addition, 
we see that increased luminosity at 8 TeV is unlikely to significantly improve the reach 
of the "vanilla" searches. The impact of these null searches on the SUSY sparticle 
spectrum is discussed in detail and the implications of these results for models with 
low fine-tuning, a future lepton collider and dark matter searches are examined. 
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1 Introduction 



The LHC was designed to explore electroweak symmetry breaking and search for physics 
beyond the Standard Model (SM). The recent discovery of the Higgs boson [l] fulfills one 
of these goals. This discovery is a triumph for science as it profoundly deepens our under- 
standing of the universe and completes our picture of the SM. However, many fundamental 
questions remain outstanding, and it is imperative that the LHC experiments continue to 
search for new physics that may shed some insight or, possibly, open new doors of investi- 
gation. 

To date, the onslaught of LHC data has challenged our most cherished theories of 
new physics [2] . Search after search further constrains the parameter space of many models, 
with some scenarios now being excluded. In addition, the property of naturalness is now 
called into question. Either new physics appears at the electroweak scale, or Nature has 
tuned the mass of the Higgs boson to one part in ~10^^. 

Supersymmetry (SUSY) elegantly solves this hierarchy problem and is the most mo- 
tivated and well-studied theory beyond the SM [s]. Naturalness considerations allow us to 
set limits on the values of parameters, given an upper bound on the level of fine-tuning 
permitted. In our recent work |4], we showed that employing the Barbieri-Giudice defini- 
tion fs'ilG] of fine-tuning suggests the existence of at least one light stop (below ~1 TeV), and 
implies that the higgsino mass and stop mixing l^^l are below ~456 GeV and ~2.2 TeV, 
respectively. On the other hand, a Higgs mass of ~126 GeV is difficult to achieve in the 
Mininal Supersymmetric Standard Model (MSSM) and requires a large stop mass or large 
stop mixing to drive large radiative corrections. This spells trouble for the simplest or con- 
strained SUSY models, such as mSUGRA, which have difficulty satisfying all experimental 
constraints (including the Higgs mass) while meeting the requirements for low fine-tuning. 
In particular, ATLAS and CMS have placed bounds (7||8] of ~1.2 TeV on the gluino and 
first and second generation squark masses in these classes of models, with the result that 
typical measures of fine-tuning in these scenarios are soaring to the level of one part in 10^ or 
higher. This has led to the claim that the constrained MSSM (cMSSM) paradigm is "hardly 
tenable" lol at this time. 



Given this situation, it is imperative to explore the parameter space of supersymme- 
try in more detail, in order to leave no corner unturned. The p(henomenological)MSSM is 
particularly well-suited to this task. In this scenario, no specific theoretical prejudice is in- 
troduced at the GUT scale, or associated with a SUSY breaking mechanism, and a short list 
of experimentally motivated considerations reduce the MSSM to 19 (or 20) real, weak-scale 
parameters, corresponding to a neutralino (or gravitino) Lightest Supersymmetric Particle 
(LSP) [To]. This allows for the study of LHC results in a more general fashion, and con- 
siders mass patterns and signatures which are not possible in the cMSSM and other more 
conventional SUSY scenarios. The pMSSM has thus garnered much attention in the liter- 
ature [4, 11 13 . In particular, we have recently shown [4,12 that the 7 TeV LHC SUSY 
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searches allow for pMSSM models with light gluinos and squarks, and that a 126 GeV Higgs 
is easily accommodated. In addition, we found a handful of models in a specific corner of 
parameter space with light stops/sbottoms and Higgsinos that allows for 1% fine-tuning or 
better. 

However, the LHC experiments have recently released a plethora of new search results 
which encompass 8 TeV updates of 7 TeV results, new signatures and new techniques. 
Here, we re-examine the pMSSM in light of this new data, and determine if our previous 
encouraging results still stack up. In particular, we investigate whether the recently released 
3rd generation LHC SUSY searches, involving top and bottom final states, exclude our 
pMSSM models with light stops/sbottoms and low fine-tuning. This provides the motivation 
for the present work, where we examine all of the ATLAS SUSY related search results released 
before mid-September 2012; this corresponds to 22 analyses with ~5 fb~^ at 7 TeV and ~6 
fb-i at 8 TeV. 



The set of pMSSM models we consider were previously generated [12] by imposing 
the following set of minimal assumptions on the R-parity conserving MSSM: (i) CP con- 
servation, [a) minimal flavor violation at the electroweak scale, {iii) degeneracy of the flrst 
and second generations of sfermion masses, [iv) the first two generations have negligible 
Yukawa couplings and A-terms. In addition it is assumed that the LSP is the neutralinc|^ 
We perform a random scan employing fiat priors over the following parameter ranges: 100 
GeV < mL,3,ei3 < 4 TeV, 400 GeV < mQ,,^,,^, < 4 TeV, 200 GeV < mQ,,u,,d3 < 4 TeV, 
50 GeV < |Mi| < 4 TeV, 100 GeV < |M2,/i| < 4 TeV, 400 GeV < IM3I < 4 TeV, \Ar,t,b\ < 
4 TeV, 100 GeV < < 4 TeV, and 1 < tan/3 < 60. We generated 3 million points 
and subjected them to the global data set, including, (i) theoretical considerations of no 
tachyons or color/charge breaking minima as well as stable vacua, (ii) precision electroweak 
data (p-parameter, invisible width of the Z and the W boson mass), {Hi) heavy fiavor physics 
(6 —7- 57, B —7- Tu, and Bs — )■ yU/i), (iv) collider searches (LEP and Tevatron direct Higgs and 
SUSY searches, and LHC stable particle searches), and (f ) astroparticle physics constraints 
(dark matter direct detection and the WMAP relic density, which we employ only as an 
upper bound). Roughly 225k models survived these constraints and form our core pMSSM 
model sample [l2]. Of these, approximately 45k models contain a lightest Higgs boson mass 
in the range 126 ± 3 GeV and form what we will refer to as our "Higgs subset" For 
further details, please see [4,12 



In what follows, we will find that even with more energy and more searches at the 
LHC, the pMSSM retains the favorable properties listed above. Specifically, we will see that 
light (< 600 — 700 GeV) 1st /2nd generation squarks and gluinos are consistent with the data, 
light (< 400 GeV) stops/sbottoms can be present, and most of our low fine-tuning models 
remain viable. In addition we will see that increased luminosity at 8 TeV results in only 
marginal improvements in the reach of the standard "vanilla" searches. We will examine 



■'■We have also generated a large set of pMSSM models -with a gravitino LSP, but these have some distinct 
collider signatures and ■will be considered elsewhere. 
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the causes and consequences of these results. In particular, we will see that light Higgsinos 
can be kinematically accessible at the ILC, and that the next generation of dark matter 
direct detection experiments have sensitivity to the low fine-tuning parameter region. In the 
next section we describe the recent LHC SUSY analyses that we incorporate, as well as the 
necessary modifications we implemented in the standard fast detection simulation PGS |14j . 
In Section 3 we present our results and their implications before concluding. 



2 LHC SUSY Searches 

In this section we describe the LHC search channels employed in our study and the corre- 
sponding implications of our fast detector simulations. 



2.1 Vanilla, Multilepton, and 3rd Generation Searches 

We begin by presenting the LHC analyses considered for this study, where we closely follow 
the ATLAS search strategies. Generally, we have included every ATLAS SUSY search publi- 
cally available as of September 2012 that is expected to significantly constrain the neutralino 
LSP pMSSM model set described above. An exception is the searches involving final-state 
taus, due to the current limitations of the fast detector simulator PGS jl4] which we em- 
ployed. In order to apply the ATLAS analyses, we have made numerous modifications to 
the PGS code, aiming to mimic the experimental environment as closely as possible. Details 
on our general strategy, and some of these changes, may be found in jl2]. Here, we describe 
the additional modifications that were necessary to reproduce the ATLAS searches for this 
work. 



First, we examine the ATLAS ~5 fb ^ 7 TeV jets plus MET p^, multi-jets fl6|, and 



one lepton |17 searches. These were considered previously [12j, but remain relevant even 
after the 8 TeV data, as will be seen below. We will refer to these searches as the "vanilla" 7 
TeV SUSY analyses, as opposed to the more specialized search channels, with signal regions 
involving third generation quarks or additional leptons. 

Now that first and second generation squarks have been excluded at low mass in 
most standard scenarios, there has been much recent interest in searching explicitly for light 
third generation squarks. Such squarks would, of course, be particularly appealing for a 
supersymmetric resolution of the hierarchy problem. Several LHC analyses have sought 
evidence of direct stop or sbottom production by considering final states with 6-quarks, both 
with and without leptons, and we have implemented many of them in this work. The various 
analyses target different squark mass ranges. As our scan ranges for the third generation 
squark mass parameters start at 200 GeV, we have not considered searches that are directed 
towards stops or sbottoms lighter than the top quark mass. Examining one such search |18|, 
we have verified that it is very unlikely that a pMSSM model with a third generation squark 
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above the top quark mass would be excluded, yet missed by more conventional searches 
for heavier third generation squarks. With this in mind, we have implemented all relevant 



ATLAS searches for third generation squarks heavier than the top quark 19 



The third generation searches use a variety of 6-tagging algorithms at different ef- 
ficiency points, reflecting both the rapid evolution of 6-tagging methods and the higher 
6-tagging efficiencies required by recent searches involving multiple 6-quarks. We have re- 
moved the original PGS 6-tagging efficiency functions, which are based on the CDF Run II 
detector, and instead employ the built-in PGS "truth" jet type function for each jet, which 
outputs the type of the highest energy parton within 20 degrees of a jet. When testing 
whether an event passes the cuts for a given signal region, we take the true type of each 
jet in the event and apply an appropriate 6-tagging efficiency function or light quark/gluon 
rejection factor. The efficiency functions and rejection factors vary for different signal regions 
and on the 6-tagging algorithm chosen by ATLAS; we estimate these from fits to ATLAS 
data [20|. As an example. Figure [l] shows the efficiency function we applied to true 6-jets 
in order to reproduce the MVl algorithm at an operating point that gave 60% efficiency 
on a tt sample. By fitting the ATLAS data with a logarithmic function for each 6-tagging 
algorithm at each operating efficiency used, we are able to match the experimental 6-tagging 
procedure well. 

There are also LHC searches using 7 TeV data with signal regions involving multiple 
leptons. These are directed at both colored sparticle production with subsequent decays 
through electroweak gauginos, as well as direct slepton/gaugino production. In principle, 
both of these channels are possible in the pMSSM, where the sparticle spectrum has no a 
priori structure. Again, we do not include searches involving taus because PGS does not 
currently realistically reproduce tau efficiencies and fake rates at the LHC. However, all 



current ATLAS analyses involving electrons and muons are implemented 21 



In addition, searches have already been reported using the 8 TeV LHC data from 
2012. In particular, ATLAS has performed analyses for final states with jets plus MET [s]. 



many jets 22 , one lepton with jets plus MET 23 , and same-sign dileptons with jets plus 
MET 24 , using ~6 fb~^ of 8 TeV data. Many of the signal regions in these analyses are 
quite similar to those of the corresponding 7 TeV searches, though some cuts are slightly 
tightened. As we will see in Section |3.1[ these small changes for the 8 TeV signal regions 



have the effect of reducing sensitivity to some models relative to the older 7 TeV searches. 



As in 12 , we also extrapolate the 8 TeV analyses to 25 fb~^ of integrated luminosity, asking 
how many models would be excluded if the searches were repeated with the same cuts but 
more data. For this extrapolation, we scale the background and corresponding error in each 



signal region and calculate the signal event limit using the CL^ method 25 , assuming that 
the observed number of events at 25 fb~^ matches the expected background. Because of 
downward fluctuations in the observed numbers of events in many 8 TeV signal regions at 6 
fb~^, the estimated sensitivity with 25 fb~^ can be comparable to, or even worse than, that 
obtained at 6 fb~^. In fact, we will find that the fraction of models in our set that would be 
expected to be excluded with 25 fb~^ of data at 8 TeV is not significantly higher than the 
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Figure 1: 6-tagging efficiency as a function of pt for the MVl algorithm. The points show 
ATLAS data obtained using a ti sample at an operating point where the algorithm's overall 
efficiency was 60%. The curve is our fit to the ATLAS data. We use the fit function on true 
6-jets to implement 6-tagging for signal regions using this algorithm and operating point. 
Similarly good ffis are obtained for other algorithms and operating points. 
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fraction of models that is already excluded. 



2.2 Long-Lived Sparticles 

In addition to the standard searches using jets, leptons, and missing energy, searches for 
more exotic detector signals, particularly those associated with long-lived charged particles, 
provide important limits on supersymmetric models. For models with neutralino LSPs, the 
lightest chargino is frequently metastable when the LSP is a nearly pure wino eigenstate. 
In this case the next-to-LSP (NLSP) is usually a wino-like chargino, and the mass splitting 
between the two is typically only ~160 MeV, leading to detector-scale decay lengths. The 



prevalence of nearly pure wino LSPs in our model set 12 results in a significant fraction 
(~4%) of models in which the lightest chargino has an unboosted decay length over 7m 
{e.g., long enough to traverse the CMS muon spectrometer), and even more models (~25%) 
in which the lightest chargino decays in the tracker, calorimeter, or muon chamber. In a 
few cases, sparticles other than the xf cire metastable as a result of small mass splittings 
with the LSP. However, these cases are quite rare since they rely on accidental degeneracies 
that are uncommon due to the random nature of our parameter scan. As a result, only 
a tiny fraction (0.14%) of our models have charged sparticles that are not charginos and 
have unboosted decay lengths beyond 7m. Our analysis is thus tailored to chargino NLSPs, 
although it is applicable to any uncolored metastable sparticle. Models with metatable 
colored sparticles require modifications to hadronization and decay routines within PYTHIA 
and will be considered in our study of the collider signatures of the pMSSM with a gravitino 
LSP |26|. 

Metastable charginos produce several different signatures depending on the location 
of their decay vertex. Charginos decaying inside the tracker produce an extremely soft 



pion, resulting in a disappearing track which is the focus of a dedicated ATLAS search 27 
Charginos decaying outside of the tracker will be detected by searches for highly ionizing 
particles, using either the inner detector alone (if the chargino decays before leaving the 
muon spectrometer) or both the inner detector and muon chamber. Both ATLAS and CMS 
search for highly ionizing particles; we choose to apply the CMS search j28| to our model set 
as the variety of signal channels and efficiencies presented are helpful for validation of our 
simulation. 

Metastable charginos require careful treatment in several stages of our analysis. As 



described in 11 , we employ analytic formulae to calculate chargino decays for models in 
which the lightest chargino and LSP masses are split by less than 1 GeV. We have also found 
it necessary to modify the PCS package. Within PCS, we have altered, or added, func- 
tions reconstructing a metastable charged particles's momentum, velocity, and calorimeter 
response, and have modified the output to include the location of any decay vertex, the pres- 
ence or absence of a muon spectrometer track, the tracker and calorimeter isolation sums, 
and the particle velocity. Our analysis routine uses the particle velocity to compute the 
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ionization energy loss in the tracker before applying cuts specific to the individual search. 
These modifications are described in detail below. 



The most sensitive searches for very long-lived charginos rely heavily on the muon 
system. We therefore introduce several modifications to the standard PGS treatment of 
the muon spectrometer response to a heavy metastable charged particle. First, we have 
implemented a special momentum resolution function for particles reconstructed in the muon 
spectrometer with > 100 GeV. This is necessary because the standard spectrometer 
momentum resolution function in PGS is simulated by the momentum resolution of the 
inner tracker, which is outperformed by the muon spectrometer at high p^. Our muon 
spectrometer momentum resolution function is a copy of that used to simulate the inner 
detector momentum resolution, but employs input parameters mimicking the ATLAS muon 
spectrometer. (The difference in resolution between the ATLAS and CMS muon systems 
is not important for our purposes.) The main motivation for including a dedicated muon 
resolution function is to decrease the extraneous MET resulting from mismeasurement of 
long-lived sparticles that frequently have a very high p^. Second, we have strengthened the 
requirements for track reconstruction within the muon spectrometer (which was previously 
only a pseudorapidity cut). A metastable charged particle can fail reconstruction in the 
muon spectrometer by decaying before it produces enough hits for reconstruction. Since the 
requirements for reconstruction depend on the detailed geometry of the muon spectrometer, 
we conservatively require muon candidates to decay outside of the muon spectrometer (r > 
7 m in the rapidity region of interest). Additionally, slow-moving particles may produce 
delayed hits in the muon spectrometer that are not associated with the primary vertex. We 



therefore subject the track reconstruction to the velocity-dependent efficiency given in 29 



fa 



for /3 < 0.3, 
2.5/3 - 0.75 for 0.3 < /3 < 0.7, 

1 for /3 >0.7, 



where face is the fraction of stable charged particles with a given velocity and a rapidity within 
the muon spectrometer acceptance that are reconstructed by the muon system. Finally, 
timing measurements within the muon spectrometer provide a measurement of /3~^ by using 
time-of-flight. We simulate this measurement by smearing the truth-level value of with 
a Gaussian with a width of 0.06, as given by [28]. 

Although muons and muon-like particles typically deposit only a few GeV of energy 
in the calorimeters, very slow-moving particles (/3 < 0.1) can deposit a significantly larger 
amount of energy. Since few metastable charginos are produced with such a low velocity, 
however, we simply neglect the calorimeters and use the muon track to correct the missing 
energy measurement. While this approach slightly overestimates the amount of missing 
energy in events with slow-moving metastable particles which fail to produce a muon track, 
the effect on our results is insignificant. 
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The CMS search for heavy stable charged particles uses a modified version of the 
Smirnov-Cramer-von Mises discriminant, denoted by las, to estimate the compatibility of 
the observed charge deposition for a given track and that expected for a particle near the 
minimum of ionization. The discriminant is defined by the following formula: 

i=l 

where J is the number of hits and Pi is the sorted probability for the energy deposited by 
a minimum-ionizing particle to be less than or equal to the energy deposited in the i*^ hit. 
Our analysis code uses our modified PCS output to calculate las in the following manner. 
First, the most probable ionization is calculated using an experimental parametrization: 



with K = 2.559 and C = 2.772 



MeV 



28 



This expression agrees closely with the 



Bethe-Bloch formula in the range 0.4 < l3 < 0.9; for values of j3 below (above) this range, 
the precise details can be neglected because the particle is clearly highly ionizing (minimum 
ionizing). Second, we assume that the particle has 16 hits in the CMS tracker l30], and choose 
an energy deposition (^)h,j for each hit from a Gaussian of width 0.559 centered at 
(^)^p. (The width is determined by a Gaussian fit to the ionization energy distribution 
for minimum bias events in Figure 4 of (STj.) We model the saturation of the CMS pixel 
detector with a threshold of 9.436 by assigning the threshold value to any hit above 
threshold. Third, for each hit we calculate the probability that a minimum-ionizing particle 
would deposit a smaller or equivalent amount of energy by integrating the probability density 
function (PDF) for its ionization energy deposition from zero to (^)hit^ Finally, we sort 
the individual probabilities and sum them according to the expression given above. 

The CMS search for heavy stable charged particles (HSCPs) [28] requires an inner 
detector track with 11 hits in the silicon tracker. We approximate this high number of re- 
quired hits by simply requiring a HSCP candidate to decay outside of the silicon tracker, r 



> 1.1 m 32 . We simulate the preselection process by imposing cuts on MET, candidate 
pseudorapidity, and calorimeter and tracker isolation sums. For the tracker+time-of-fiight 
analysis, we additionally require the presence of a track in the muon spectrometer. The 
final selection is divided into search regions, defined by cuts on the values of p^,, las, the 
reconstructed mass, and (for the tracker+time-of-fiight analysis) /3. We simulate the recon- 
structed mass cut by applying a given search region only to candidates with a truth-level 



^We take the PDF from Figure 4 of 31 . While this distribution is generated from pixel detector hits for 
a sample of minimum bias events at y/s = 0.9 TeV, most of these hits come from pions near the minimum 
of ionization. As a cross-check, the PDF for minimum ionizing particles determined in a test run |33j is 
similar to the one we use. We also expect that the PDFs for the silicon strip tracker and pixel detector are 
equivalent within our desired accuracy. 
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Benchmark 


a (us) 


A (us) 


a (CMS) 


A (CMS) 


100 GeV fi 


41 fb 


20% 


38 fb 


11% 


308 GeV n 


0.42 fb 


59% 


0.35 fb 


39% 



Table 1: Validation of our simulation of the CMS tracker-only HSCP search, comparing our 
results for the production cross section and acceptance with that of CMS. 



mass at least 1 standard deviation above the cut value. This implies that an insignificant 
fraction of signal events will be removed by the mass cut. We note that all of our stable 
particles pass the mass cut for at least one search region, and even low-mass search regions 
have very high sensitivity. We validated this analysis by comparing our simulation with two 
benchmarks: Pair-produced stable staus with masses of 100 GeV and 308 GeV. The results 
of our validation are shown in Tables [T] and [2j Comparing our results with those of CMS 
shows that we reproduce the production cross section for both benchmarks. The acceptance 
of the tracker -|-time-of-flight analysis is also reproduced quite nicely, while our acceptance for 
the tracker-only analysis (particularly for the 100 GeV stau) is significantly higher than that 
predicted by CMS. Since the tracker-only analysis is only relevant for charginos decaying 
in the calorimeter or muon spectrometer, the tracker-|-time-of-fiight analysis is considerably 
more important for our model set. We note that the discrepancy seen for the light stau 
model is most likely because the light staus tend to have /5 ~ 1, making them difficult to 
distinguish from minimum-ionizing particles and introducing a strong sensitivity to the de- 
tails of ionization energy deposition. Fortunately, models with very light stable charginos 
are excluded by the LHC, the Tevatron and LEP (and therefore are excluded during model 
generation). The accuracy of our treatment is significantly better where it is most relevant, 
namely for heavier charginos with masses near the the exclusion boundary. 



The ATLAS disappearing track search 27 looks for charginos which are reconstructed 



in the silicon tracker (SCT), but not in the outer transition radiation tracker (TRT). We 
therefore require candidates to decay between the outer edge of the SCT (520 mm) and the 
inner edge of the outer TRT (863 mm). The high number of hits required in the SCT and low 
number required in the outer TRT support these sharp cutoffs. We approximate the track 
isolation requirement (no tracks with > 0.4 GeV within AR < 0.1) by using AR < 0.25; 
a negligible difference between the two requirements was observed during comparison with 
benchmarks. We then require candidates to be the highest isolated track in the event 
and impose cuts on the standard physics objects (leptons, jets, and MET). We validated our 
analysis with the ATLAS benchmarks (100 and 200 GeV charginos with 1 ns lifetimes); the 
results are presented in Table [3j Our pair-production cross sections and signal efficiencies 
for both models are in agreement with the ATLAS values, indicating that our geometric 
acceptance cuts and isolation requirements accurately reproduce the ATLAS analysis. 
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Benchmark 


a (us) 


A (us) 


a (CMS) 


A (CMS) 


100 GeV fi 


41 fb 


21% 


38 fb 


19% 


308 GeV n 


0.42 fb 


63% 


0.35 fb 


55% 



Table 2: Validation of our simulation of the CMS tracker+time-of-flight HSCP search, com- 
paring our results for the production cross section and acceptance with that of CMS. 



Benchmark 


a (us) 


A (us) 


a (ATLAS) 


A (ATLAS) 


100 GeV Xi 
200 GeV Xi 


13100 fb 
819 fb 


0.070 % 
0.11 % 


14400 fb 
808 fb 


0.066 % 
0.12 % 



Table 3: Validation of our simulation of the ATLAS disappearing track search, comparing 
our results for the production cross section and acceptance with that of ATLAS. 

3 Results and Implications 

In this section we will discuss the results of the various SUSY searches and their implications 
for the parameter space of our pMSSM model set with a neutralino LSP. 



3.1 Search Results 

The first question we address is how well these searches perform, both individually and when 
combined, in covering the set of pMSSM models that we generated. We begin by considering 
the "vanilla," generalized jets plus MET, searches at 7 and 8 TeV. Tables |4] and [5] show the 
fraction of our model set excluded by these two suites of searches, individually and when 
they are combined. Note that results are shown for both the full model set and for the 
subset of models satisfying m/^ = 126 ± 3 GeV, henceforth referred to as the "Higgs subset" . 
There are several things to observe about these results: (i) Individual models are more than 
likely to be excluded by more than one of these searches as is clear when we compare the 
combined results to those for the individual searches themselves, {ii) The 2-6 jets plus MET 
analysis is by far dominant and provides most of the model coverage at both center of mass 
energies. This is different from the ATLAS results for their corresponding cMSSM SUSY 
searches §|T5}jT7)|22}]24] ; there the coverage provided by the 1-lepton plus jets search channel 



was found to be relatively comparable to that of the 2-6 jets plus MET channel 34 . There 



are several reasons why the leptonic searches are degraded in the pMSSM case, the most 



important being the lack of light binos in the spectrum. As we have discussed earlier 12 
most of the LSPs in our pMSSM model set are either nearly pure wino or Higgsino due to 
DM constraints and so tend to have 2 or 3 nearly degenerate states at the bottom of their 
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Search 


Reference 


Full Model Set 


Higgs Subset 


2-6 jets 


ATLAS-CONF-2012-033 


21.04% 


18.53% 


multijets 


ATLAS-CONF-2012-037 


1.61% 


1.34% 


1-lepton 


ATLAS-CONF-2012-041 


3.16% 


2.80% 


Total 




21.19% 


18.64% 



Table 4: Fraction of our pMSSM models with neutralino LSPs excluded (in per cent) by 
the general "vanilla" MET ATLAS searches at the 7 TeV LHC with 4.7 fb~^ of integrated 
luminosity for both the full model set as well as for the "Higgs Subset" satisfying the Higgs 
mass constraint, nih = 126 ± 3 GeV. 



spectra. On the other hand, in the cMSSM the LSP is quite commonly bino-like (with no 
other states very close in mass) with a somewhat heavier wino. This then allows for wino 
decays through gauge bosons which subsequently produce leptons with sufficient p^, to pass 
analysis cuts in the cMSSM, whereas the corresponding leptons in the pMSSM will be too 
soft. In the subset of pMSSM models with relatively light bino-like neutralinos we indeed find 
an improved efficiency for model coverage by the leptonic channels. We might also expect the 
performance of leptonic channels to improve significantly with increasing luminosity, which 
will increase the sensitivity of these channels to final states with low branching fractions, 
e.g. leptons produced in cascade decays. (Hi) We note that employing the Higgs mass 
constraint does not significantly alter the sensitivity of the search channels in a qualitative 
manner, although there is some degradation (~ 10%) observed in the coverage provided by 
each of the searches individually as well as in the combination, (iv) Comparing the results 
obtained at 7 TeV with ~5 fb~^ and those at 8 TeV with ~6 fb~^, we see that the 8 TeV 
data increases the total amount of model coverage by ~25% and that the model coverage 
provided by the searches individually improves. This implies that more models are being 
excluded by multiple searches at 8 TeV than are at 7 TeV. 

Interestingly, a small number of models (~1.6k) are found to be excluded by the 7 TeV 
searches but not by the corresponding channels at 8 TeV. This is likely due to the increased 
jet and MET requirements of the 8 TeV searches that are now allowing models with 
somewhat degenerate spectra to be missed. There is a valuable lesson to be drawn from this 
result. As the LHC energy is further increased, it will be likely that some reasonable fraction 
of pMSSM model points with spectrum degeneracies will be bypassed by the analyses due 
to the continual strengthening of analysis cuts. For more complete model coverage it is then 
clearly useful and necessary to combine analyses that are performed at different center of 



mass energies, or to perform targeted searches for compressed spectra 35 . In our case of 
interest, however, combining the 7 and 8 TeV generalized analyses we find that the total 
model coverage increases by only ~1%. 
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Search 


Reference 


Full Model Set 


Higgs Subset 


2-6 jets 


ATLAS-CONF-2012-109 


26.51% 


23.82% 


multijets 


ATLAS-CONF-2012-103 


3.31% 


2.84% 


1-lepton 


ATLAS-CONF-2012-104 


3.30% 


3.07% 


SS dileptons 


ATLAS-CONF-2012-105 


4.88% 


4.50% 


Total 




26.90% 


24.16% 



Table 5: Same as Table |4] but now for the 8 TeV, 5.8 fb ^ ATLAS searches. 



Given these results, we can ask how well these generalized jet plus MET searches at 8 
TeV will perform when an expected luminosity of 25 fb~^ is reached by the end of the 2012 
run. To do this kind of extrapolation we need to make several assumptions: (a) the ATLAS 
analyses are unchanged from their present form in terms of cuts, etc. (although these are 
likely to be strengthened); (b) we can simply scale the SM backgrounds found by ATLAS 
in the lower luminosity data to the case of 25 fb~^ and (c) the number of events observed 
in each analysis channel is exactly matched to the corresponding anticipated backgroundsj^] 
With this set of conservative assumptions we obtain the results shown in Table [6j Here we 
see that whereas for the two leptonic searches pMSSM model coverages are increased relative 
to the lower luminosity results, those for the 2-6 jets plus MET search are actually degraded 
whereas those for the multijet final state are essentially unchanged. The reason for this is the 
result of our assumption (c) above and the fact that ATLAS saw fewer events than expected 
in many of their 2-6 jets plus MET search regions with 5.8 fb~^. Clearly the exact coverage 
will be quite sensitive to the actual observed signal vs. background ratio and this situation 
will need revisiting subsequent to the publication of the complete ATLAS 25 fb~^ results 
sometime next year. Overall, given our assumptions, it seems that higher luminosities may 
not cover too much more of the parameter space. What is really needed for more coverage 
is to go to higher energies, ^/s = 13 — 14 TeV. 



We now turn our attention to the large set of Heavy Flavor (HF) and Multilepton 
(ML) searches which, to date, have only been reported by ATLAS for the ^/s = 7 TeV 
data set, mostly with 4.7 fb~^ of integrated luminosity. As noted above, these searches 
are especially important as we would expect natural SUSY models to have relatively light 
stops/sbottoms as well as some light Higgsino-like (and possibly wino-like) gauginos. To 
date, 11 of these HF/ML searches, as shown in Table [7| have been reported by ATLAS. 
As noted above, in order to analyze these HF channels, it was necessary to modify the 
b-tagging routines employed by PGS [li] to more accurately reflect those of the ATLAS 
20 . The results shown in Table ItI display the pMSSM coverage for each separate 



detector 



^If ATLAS is "lucky" and observes fewer events than anticipated (as frequently happens) then stronger 
constraints will naturally result from the relevant analysis. 
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Search 


Reference 


Full Model Set 


Higgs Subset 


2-6 jets 


ATLAS-CONF-2012-109 


25.27% 


22.68% 


multijets 


ATLAS-CONF-2012-103 


3.31% 


2.84% 


1-lepton 


ATLAS-CONF-2012-104 


3.83% 


3.57% 


SS dileptons 


ATLAS-CONF-2012-105 


7.45% 


6.96% 


Total 




26.05% 


23.49% 



Table 6: Same as Table |5jbut now extrapolated to 25 fb ^ of integrated luminosity at 8 TeV. 

analysis, as well as for the two combinations of the set of HF and ML analyses. Here, 
again, the results are shown for both the full model set as well as for the subset of models 
with the Higgs mass constraint imposed. As in the case of the generalized MET searches, 
we see that applying the Higgs mass cut causes some small degradation (~10-15%) in the 
sensitivity of these analyses as well as in the combinations. Clearly some of these searches 
will be more powerful at probing the pMSSM than are others. It is no surprise that the 
coverage is rather poor for the very light stop search (since all our stops are required to 
be more massive than the top quark in the model generation process), as well as for the 
direct stop in the natural GMSB search. Also, since the cross sections for direct gaugino 
production are rather small and, in most of our models, these states are quite degenerate 
due to their being relatively pure wino/Higgsino eigenstates, the direct gaugino search is 
also generally not very effective at probing our model set. Other searches, such as those 
for heavy stops, direct sbottoms, and 1-2 leptons plus jets, perform significantly better and 
provide reasonable model coverage. Since these searches are relatively orthogonal to the 
generalized jet plus MET searches discussed above, we would expect that HF/ML searches 
will add significantly to the overall pMSSM model coverage and this is indeed the case. In 
particular, we find that 0.96%, 0.08%, and 1.03% of models are excluded by heavy flavor 
searches, multilepton searches, and their combination (respectively), but are not covered by 
the 7 and 8 TeV "vanilla" searches. 

In addition to searches with MET, those looking for heavy stable charged particles 
(HSCP) or disappearing tracks (DT) will be quite important for our pMSSM model set due 
to the rather large number of long-lived charginos present in these models. Such states will 
be either directly produced or may be present in low-MET cascades. In the latter case we 
can access long-lived charginos with masses beyond those accessible via direct production. 
Again the relatively frequent appearance of long-lived charginos in our model sample is due 
to the rather pure wino- and Higgsino-like nature of most LSPs and the relative degeneracies 
of the charginos and the LSP which occur in such cases. Table [8] shows the sensitivity of 
these searches for the full model set and for the Higgs subset. Interestingly, unlike the other 
searches that were discussed above, that for HSCP performs slightly better in the subsample 
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Search 


Reference 


Full Model Set 


Higgs Subset 


Gluino — 7- Stop/Sbottom 


1207.4686 


4.92% 


4.54% 


Very Light Stop 


ATLAS-CONF-2012-059 


<0.01% 


<0.01% 


Medium Stop 


ATLAS-CONF-2012-071 


0.32% 


0.24% 


Heavy Stop (01) 


1208.1447 


3.66% 


3.15% 


Heavy Stop (11) 


1208.2590 


1.94% 


1.69% 


GMSB Direct Stop 


1204.6736 


<0.01% 


0% 


Direct Sbottom 


ATLAS-CONF-2012-106 


2.47% 


2.20% 


3 leptons 


ATLAS-CONF-2012-108 


1.05% 


0.92% 


1-2 leptons 


1208.4688 


4.11% 


3.61% 


Direct slepton/gaugino (21) 


1208.2884 


0.11% 


0.10% 


Direct gaugino (31) 


1208.3144 


0.33% 


0.27% 


HF Total 




7.26% 


6.48% 


ML Total 




4.29% 


3.78% 



Table 7: Same as Table 4 but now for the HF and ML searches. 
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Search 


Reference 


Full Model Set 


Higgs Subset 


HSCP 


1205.0272 


4.03% 


4.14% 


Dis. Tracks 


ATLAS-CONF-2012-111 


2.59% 


2.21% 




[36 


39 




2.70% 


5.58% 


A/H T+T- 


1202.4083 


0.07% 


0.03% 


AU Searches 








33.89% 


33.45% 



Table 8: Same as Table |4| but now for the non-MET searches. The corresponding combined 
results obtained from all searches is also shown. 



of models where the Higgs mass constraint is applied. 

Finally, before we combine the results from all these searches to extract the total 
pMSSM coverage at the LHC, we must recall the significant contributions arising from both 



the searches for Bs — )■ jjL^ (i 36 -39] as well as for A, H ^ t^t 40 which were discussed 



in detail in our earlier work [12] . The former, in particular, is relatively powerful, and can 
allow us to access model points with heavier gluino and/or first and second generation squark 
masses than do the present MET searches. These results are reproduced in Table |8] where 
we see that for the subset of models where the Higgs mass constraint has been enforced the 
Bs — 7- is significantly more effective. 

Combining all these various searches (from Tables |4| |5l [T] and [s]) together without 
(with) the Higgs mass constraint imposed, we see from Table |8] that 33.89 (33.45)% of our 
pMSSM models are now excluded by the combined data from the LHC. Note that when all 
the searches are combined, the overall effect of the Higgs mass constraint is very weak. To 
make a significant improvement in these numbers the complete search results for 8 TeV will 
need to be included. Of course, going to a/s ~ 14 TeV will provide a major step in model 
coverage. 



3.2 Implications 

We now turn to a discussion of some of the implications of these various searches. The first 
question we address is how the distributions of the various sparticle masses are influenced by 
the set of negative SUSY search results obtained so far at the LHC. Furthermore, we would 
also like to know if some mass ranges in the pMSSM are now entirely excluded by the LHC 
searches. These questions are answered by examining the nested set of histograms presented 
for the various sparticle types in Figures |2]j4j 
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Figure 2: Histograms showing the distribution of gluino (upper panel) and hghtest lst/2nd 
generation squark (lower panel) masses in our pMSSM model set with a neutralino LSP. The 
effect of sequentially applying LHC searches on these distributions is shown as a series of 
colored lines in the following order from top to bottom: Full Model Set (black), 7 and 8 TeV 
Vanilla Searches (red), Heavy Flavor (green), Multileptons (blue), HSCPs and Disappearing 
Tracks (magenta), Bg — )■ /^"""/i" and H/A — t- r+r~ (cyan), and m^j = 126 ± 3 GeV (brown). 
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Figure 3: Histograms showing the distribution of the hghtest stop (upper panel) and hghtest 
sbottom (lower panel) masses in our pMSSM model set with a neutralino LSP. The line 
colors are as in Figure [21 
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Figure 4: Histograms showing the distribution of the hghtest neutrahno (upper panel) and 
hghtest chargino (lower panel) masses in our pMSSM model set with a neutralino LSP. The 
line colors are as in Figure [2j 
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In these figures we see the distribution for the masses of various sparticles for models 
that remain valid after evading a series of LHC experimental searches, with each line indi- 
cating the effect of applying an additional search in the following order: 7 and 8 TeV Vanilla 
Searches, Heavy Flavor, Multileptons, HSCPs and Disappearing Tracks, Bg — )■ fi^fi~ and 
H/A T+T-, and m^ = 126 ± 3 GeV. 

Let us first consider the case of the gluino. In Figure |2] we see, as expected, that 
the 7 and 8 TeV general MET searches take the largest bite out of the gluino mass distri- 
bution, particularly at the low mass end. The other searches tend to remove model points 
corresponding to heavier gluino masses which are not directly kinematically accessible. In- 
terestingly, and we will discuss this further below, we see that a sizable fraction of models 
still exist with gluino masses below ~700 GeV. In the case of the lightest first/second gener- 
ation squarks, while the 7 and 8 TeV general MET searches remove many models from the 
distribution, they are most effective below ~900 GeV while at larger masses the constraints 
are indirect. The HSCP searches are seen to be more effective in this case particularly in the 
intermediate mass region ~1.0-1.5 TeV. We again see that models with somewhat lighter 
squarks with masses just above ~600 GeV still remain viable. 

Lighter stops/sbottoms are a common prediction in natural models with low fine- 
tuning, so it is particularly important to examine the sensitivity of the HF searches to 
these states. In Figure |3] we see the spectra for both light stops and sbottoms in the most 
interesting mass region below 1 TeV. In addition to the expected large impact of the 7 and 8 
TeV general MET searches we see, particularly in the mass range below ~ 650 GeV, that the 
HF searches are making a significant dent in the number of surviving models. However, the 
number of models with light stops or sbottoms below ~ 400 GeV still remains significant. 
When the Higgs mass constraint is included there is seen to be a serious reduction in the 
number of remaining models with lighter stops. This is not too surprising as heavy stops, 
and/or large stop mixing, is required in the pMSSM to obtain Higgs masses in the required 
126 ± 3 GeV range [4]. We note, however, that since the Higgs mass requirement does not 
directly affect the existence of lighter sbottoms, a significant number of models remain with 
relatively low sbottom masses. In fact, after applying the Higgs mass constraint the sbottoms 
are generally found to be lighter than the stops. The fact that relatively light stops/sbottoms 
remain in the spectra of the surviving models is not too surprising as the standard stop and 
sbottom searches usually assume a single decay mode dominates. As we saw in our earlier 
work |4|, specifically in the case of models with low fine-tuning, both stops and sbottoms 
have multiple decay modes with comparable branching fractions making their observation 
in any one channel significantly more difficult. This is further compounded by the multiple 
cascades that can take place due to the subsequent decays of the light electroweak gauginos 
(down to the LSP) into which the stops and sbottoms themselves decay. 

Figure |4] shows the mass distributions for both the LSP and the lightest chargino 
which are seen to track each other quite closely overall. This is to be expected as the bulk 
of our pMSSM models have LSPs which are dominantly wino- or Higgsino-like and thus 
have a nearby charged partner. The 7 and 8 TeV searches (indirectly) cause almost all of 
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the reduction in the model set and none of the other searches are seen to be particularly 
effective (except for HSCP) outside of the lowest mass bins. The effectiveness of the HSCP 
search is due to the existence of many HSCPs which are nearly LSP-degenerate charginos. 

Let us return to the models with relatively light squarks and gluinos which remain 
valid after all of the search constraints are applied. Figures [5] and [6} display how the density 
of the gluino and lightest first /second generation squark masses respond to the full set of 
searches (except for the application of the Higgs mass cut which we now know does not have 
much of an effect except for lowering the overall statistics). Figure |5] shows results for the 
entire model set distribution whereas Figure |6] highlights only the squark and gluino mass 
range below 1.5 TeV. Recall that this is essentially the excluded region in the case of the 
cMSSM due to the ATLAS jets+MET generic searches. In both Figures the top panel shows 
the distribution for the original pMSSM model set as generated, whereas the lower panel 
shows only the surviving models after all of the LHC constraints have been applied. From 
both figures it is clear that the searches have predominantly excluded models with either 
light gluinos and/or light squarks with the gluinos showing greater sensitivity as expected. 
However, the zoomed-in panels shown in Figure |6] demonstrate that models with relatively 
light squarks and gluinos do remain in the surviving model sample. In particular, models 
with gluinos below ~700 GeV or squarks below ~600 GeV remain viable. 

To understand how these models survive the onslaught of the multiple LHC searches, 
we show in Figure [7] the mass distribution histograms for the LSP in the subset of sur- 
vivor models for gluinos (top panel) or first / second generation squarks (bottom panel) with 
masses below 800 GeV as the various searches are imposed. One sees in both cases, after all 
the searches are combined (but primarily arising from just the 7 and 8 TeV general MET 
searches), that the LSP masses are quite peaked near the upper limit of their possible mass 
range. This demonstrates that the models with light squarks and/or gluinos that are sur- 
vivors are mostly the result of compressed spectra, particularly in the case of light surviving 
gluinos. There are some models with light squarks that evade the searches for other reasons. 

It is interesting to ask if the general 8 TeV searches or any of the HF or ML searches 
are able to exclude some of the 13 low fine-tuning model^ i.e., those with fine-tuning better 
than 1%, with neutralino LSPs and rrih = 126 ± 3 GeV as were discussed in our earlier 
work [4|. We find, in fact, that only one of these models (the one with the lightest stop 
and sbottom as one might expect) is actually excluded by the HF analyses, in particular, 
by both the b-jet plus MET and light stop searches. Furthermore, we find that this same 
model, as well as an additional 3 of the remaining 12 models, are excluded by the "vanilla" 
MET searches at 8 TeV as light stops and sbottoms produce jets plus MET and these models 
have (in 2 out of 3 cases) relatively light gluinos. If we consider the somewhat larger set of 
50 neutralino LSP models satisfying the more relaxed fine-tuning constraint (the fine-tuning 
parameter A < 120) with rrih = 126 ± 3 GeV, we find that 34 of them survive the 7 TeV 
general MET searches. 24 of these also survive the corresponding 8 TeV "vanilla" MET 

"Here we have employed the standard Barbieri-Giudice fine-tuning definition 5 6. 
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Figure 5: The density of models in the hghtest lst/2nd generation squark - gluino mass plane, 
before (upper panel) and after (lower panel) applying the LHC SUSY search constraints. 
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Figure 6: Close-up of the low mass region of the previous figure. 



22 



1000 



800 



600 



400 



200 




400 600 
m,« (GeV) 



1000 



3000 



2500 




1000 



Figure 7: Histograms of the lightest neutrahno mass in models with a gluino (upper panel) 
or lightest lst/2nd generation squark (lower panel) with a mass below 800 GeV. The line 
colors are as in Figure [2j 
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searches. If we now apply all of our LHC search constraints above we find that 22 models 
still survive at present. The last 2 models are excluded by the non-MET searches. It is 
interesting that the 3'''^ generation searches are not having much of an effect on this class of 
pMSSM models, likely due to the availability of multiple decay channels. 

Of course there are alternative ways of looking at fine-tuning. A recent example of 



this is provided by Refs. 41 , 42 which considered a fine-tuning measure that depends only 
upon weak scale quantities which is arguably what we are dealing with when discussing 
the pMSSM. The fine-tuning requirements in such a case are much looser than in the more 
traditional Barbieri-Giudice (B-G) measure so that many more models will have low fine- 
tuning employing this definition. Of course, the largest contributor to fine-tuning using 
either measure arises due to the size of the /i, i.e., FT^ ~ /i^/M|. In Figure [s] we show the 
electroweak fine-tuning measure presented by Baer et. al in Refs. [4T||42] for our neutralino 
pMSSM model set and for the Higgs subset. This should be compared with the results shown 
in Figure 17 of our earlier work [4] where the B-G measure is employed. As can be seen in this 
figure, the number of low fine-tuning (<100) models using this weak scale measure is several 
orders of magnitude larger, even after the Higgs mass constraint is applied. Furthermore, 
there are now a statistically reasonable subset of models with fine-tuning values <10. 

It is hopeful that a high energy e^e^ collider will be built in the future to study the 
Higgs boson and other possible new physics discovered at the LHC. It is thus important to 
ask what SUSY particles, if any, remain kinematically accessible to such a collider with a 
given value of ^/s accounting for the LHC search results. Table [o] presents the fraction of our 
pMSSM models that survive all LHC constraints (including that for the Higgs mass) which 
have a given sparticle kinematically accessible via pair production for different possible values 
of y/s for such a high energy e~^e~ collider. For ^/s = 250 GeV, which is appropriate for a 
Higgs factory, there is a small sample of models for which the first two neutralinos and the 
lightest chargino are accessible; this fraction grows substantially for ^/s = 500 GeV. These 
results should not be too surprising since most of our LSPs are either wino- or Higgsino-like. 
This suggests that light gauginos will be the most likely target of a low energy e^e~ collider 
given the present LHC data. This result is, of course, further reinforced by low fine-tuning 
considerations. In the B-G approach, the largest tree-level fine-tuning contribution arises 
from n, FT^ ~ 4/i^/M|. Thus requiring fine-tuning < 100 (10) implies that < 456 (144) 
GeV so that the light Higgsino-like states should be anticipated in the low-mass part of the 
spectrum. 

At y/s = 500 GeV, we begin to see small samples of models where both sleptons and 
the next heaviest set of electroweak gauginos become kinematically accessible, although the 
numbers are still quite small. Once ^/s = 1 TeV is reached, small samples of models with 
kinematically accessible light stops and sbottoms are present. It is clear, however, that a 
sizable fraction of the sparticle spectrum will only become accessible at collision energies 
above ~2 TeV. 

Our results also have implications for dark matter searches. For example. Figure [9] 
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Figure 8: Number of models with FT < A according to the fine-tuning measure of Baer et. 
al in Refs. (41 42 . From top to bottom the histograms are for the full neutralino LSP model 



set (green) and the neutrahno Higgs subset (blue). 
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Tabic 9: Fraction (in per cent) of our pMSSM model set with neutralino LSPs surviving all 
LHC constraints (including the Higgs mass) with kinematically accessible sparticles at an 
e+e~ hnear collider of various center of mass energies. The lack of an entry signals that the 
fraction is below 0.1%. This model sample corresponds to ~30k parameter space points. 
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shows the direct detection spin-independent cross section (scaled by the ratio of the rehc 
density to the total dark matter density) as a function of the LSP mass, and displays the 
recent XENON-100 bound [43] as well as the expected reach of XENON-IT 44 for reference. 
The upper left (black) scatter plot shows the results for our original pMSSM neutralino model 
set as generated, while the upper right panel shows the corresponding subset of models 
surviving all the LHC search constraints except for the Higgs mass requirement. The lower 
left (brown) panel includes this constraint. We see from this figure that, apart from the 
overall density of points, the LHC searches do not have much influence on the overall nature 
of the distribution of cross section predictions. The Higgs mass requirement, though thinning 
the point density significantly, is also seen to have little impact on these expectations. Thus 
we conclude that pMSSM searches at the LHC have very little impact on the possible range 
of spin-independent relic density-scaled cross sections, and that indeed these searches are 
complementary. The lower right panel shows the predicted relic density-scaled cross sections 
for the remaining low fine-tuning models discussed above. Here we see that these models 
should mostly be accessible to the XENON-IT experiment. 



4 Conclusion 

With the remarkable performance of the LHC, many searches for supersymmetry have been 
performed this year with both 7 TeV and 8 TeV data. It is important to ensure that as 
analyses are designed to search for new physics, experimental collaborations strive to obtain 
maximal coverage of the vast parameter space of SUSY. The most simple searches undoubt- 
edly leave corners of this space unexplored, necessitating the identification of such regions 
and the development of new search strategies to survey them. Here, we have continued 
our previous work in this area by studying the capabilities of the latest LHC searches to 
constrain a large set of pMSSM models. Compared to previous searches, we find that signif- 
icantly improved coverage of the pMSSM is afforded by the additional analyses. While our 
original model set did not require a Higgs mass near 126 GeV, our results are found to be 
generally independent of this additional constraint. 

The first searches for SUSY with 8 TeV data employ signal regions that are very 
similar to those in their 7 TeV counterparts. With slightly tighter cuts and several downward 
fluctuations in the numbers of observed versus expected events, these searches successfully 
exclude models that had remained viable after the corresponding 7 TeV analyses. An increase 
in integrated luminosity at 8 TeV, however, will not cause the coverage of these searches to 
improve much further, unless either downward fluctuations of events continue to occur or 
reductions in the background uncertainties can be achieved. We also caution that a small 
fraction of models that are allowed by the 8 TeV searches had produced enough signal 
events to be excluded by the analogous 7 TeV searches. This result is symptomatic of the 
unavoidable trade-off between obtaining search sensitivity in new areas of parameter space 
and retaining sensitivity in previously probed regions. 
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Figure 9: Scatter plots showing the mass of the hghtest neutrahno versus the spin- 
independent cross section for: all models (upper left), models which survive the LHC SUSY 
searches (upper right), models which survive the LHC SUSY searches and have mt = 126 ±3 
GeV (lower left), and models with low fine-tuning (lower right), which correspond to the red 
points. 
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In addition, directed searches for stops/sbottoms, sleptons, and electroweak gauginos 
involve less conventional signatures of supersymmetry, particularly 6-jets and leptons in the 
final state. While these searches rule out fewer pMSSM models overall than the more general 
SUSY searches, they enjoy success in excluding some spectra that are missed by the latter. 
Often, such models have heavier gluinos and lst/2nd generation squarks, making it difficult 
for conventional searches to see them. Because most of our gauginos are wino and Higgsino 
multiplets with very small mass splittings, the heavy flavor signal regions have comparatively 
better sensitivity to our pMSSM models than the direct gaugino searches. 

The last class of searches that we have included in this work involves signatures 
without missing energy. Given the significant population of models with stable charginos 
within our model set, it is not surprising that searches for heavy stable charged particles and 
disappearing tracks exclude significant fractions of the models. Updates to the Bg A*"*"/^" 
and A,H ^ t^t~ analyses have improved the LHC coverage of our model set as well. 
Notably, since the signal regions for these non-MET searches involve criteria that are very 
different from those of other searches, the coverage gained from these analyses is nearly 
orthogonal to that obtained from the MET-based searches. 

After including all of these analyses, we find that many pMSSM spectra with light 
sparticles survive. While supersymmetry continues to be limited by the lack of collider signa- 
tures thus far, the state of the full MSSM is far better than suggested by limits on simplified 
models or more constrained scenarios such as mSUGRA. In particular, we find surviving 
models containing lst/2nd generation squarks, gluinos, and 3rd generation squarks with 
masses below 600, 700, and 400 GeV, respectively. Models with light colored sparticles that 
are not excluded by LHC searches generally have compressed spectra. We also observe viable 
models with fine-tuning near 1%, indicating that the MSSM still has significant potential to 
provide a natural solution to the hierarchy problem. 

Finally, we have considered the potential of experiments outside the LHC to observe 
supersymmetry. A linear collider would be able to produce many of the light gaugino states 
in our viable models. We emphasize that if the MSSM is indeed natural, light Higgsinos 
will be copiously produced at a future e+e^ collider. Astrophysical experiments provide 
a complementary method of searching for supersymmetry by seeking signatures of dark 
matter, which is composed of neutralinos in the models considered here. By searching 
for direct interactions of the LSP with matter, an upgraded XENON experiment would 
probe a significant fraction of our pMSSM model set. Both of these avenues of exploration 
are orthogonal to the LHC supersymmetry searches, and would be particularly effective in 
probing natural SUSY. 

Having surveyed the vast landscape of minimal supersymmetry with the pMSSM, 
we declare it to be alive and well. We look forward to upcoming collider and astrophysical 
results, in the hope of finally revealing the hitherto elusive face of SUSY. 
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